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Abstract 

We investigate second order quasilinear equations of the form 

where u is a function of n independent variables x%, ...,x n , and the coefficients depend 
on the first order derivatives p 1 = u Xl , ■~,p n = u Xn only. We demonstrate that the natural 
equivalence group of the problem is isomorphic to SL(n + 1,-R), which acts by projective 
transformations on the space P n with coordinates p , ■■■,p n . The coefficient matrix /y de- 
fines on P n a conformal structure fij(p)dp l dp : > . In this paper we concentrate on the case 
n = 3, although some results hold in any dimension. The necessary and sufficient conditions 
for the integrability of such equations by the method of hydrodynamic reductions are de- 
rived. These conditions constitute an over-determined system of PDEs for the coefficients 
fij, which is in involution. We prove that the moduli space of integrable equations is 20- 
dimensional. Based on these results, we show that any equation satisfying the integrability 
conditions is necessarily conservative, and possesses a dispersionless Lax pair. Reformulated 
in differential-geometric terms, the integrability conditions imply that the conformal struc- 
ture fij(p)dp l dp : > is conformally flat, and possesses an infinity of 3-conjugate null coordinate 
systems parametrized by three arbitrary functions of one variable. Integrable equations 
provide an abundance of explicit examples of such conformal structures parametrized by 
elementary functions, elliptic functions and modular forms. 
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1 Introduction 



The main object of our study are second order quasilinear equations of the form 



fllU xx + f&Uyy + f^Utt + 2fl2Uxy + 2/l3« a;t + 2/ 23 % t = 0, 



(1) 



where u is a function of three independent variables x, y, t, and the coefficients depend 
on the first order derivatives u x ,u y ,ut only. Equations of this type arise in a wide range of 
applications in mechanics, general relativity, differential geometry and the theory of integrablc 
systems. Among the most familiar examples one should mention the Boyer-Finley equation, 



which is descriptive of a class of self-dual 4-manifolds [TO] , as well as the dispersionless Kadomtsev- 
Petviashvili (dKP) equation, 



also known as the Khokhlov-Zabolotskaya equation, which arises in non-linear acoustics |45] and 
the theory of Einstein- Weyl structures [2] . 

The integrability aspects of equations of the form (pQ) have been investigated by a whole 
variety of modern techniques including symmetry analysis, differential- geometric and algebro- 
geometric methods, dispersionless (9-dressing, factorization techniques, Virasoro constraints, hy- 
drodynamic reductions, etc. However, until recently there was no intrinsic approach which would 
allow a unified treatment of these and other examples. Moreover, there was no satisfactory defi- 
nition of the integrability which would (a) be algorithmically verifiable, (b) lead to classification 
results and (c) provide a scheme for the construction of exact solutions. We emphasize that 
equations of the form (pQ) differ essentially from their solitonic counterparts, and require an 
alternative approach. Such approach, based on the method of hydrodynamic reductions and, 
primarily, on the work [251 [26], see also pTJ EU E21 [TOJ [TOJ [37] , etc, was proposed in [TTJ, [20] . 
It was suggested to define the integrability of a multi-dimensional dispersionless system by re- 
quiring the existence of 'sufficiently many' hydrodynamic reductions which provide multi-phase 
solutions playing a role similar to that of algebro-geometric solutions of soliton equations. 

In Sect. 2 we briefly review the method of hydrodynamic reductions, and apply it to equations 
of the form (TjQ). This leads to a system of differential constraints for the coefficients fij, which 
are necessary and sufficient for the integrability. We demonstrate that the system of constraints 
is in involution, and prove our first main result (Theorem 1 of Sect. 2): 

• The moduli space of integrable equations of the form (TjQ) is 20-dimensional. 

In Sect. 3 we point out that the class of equations (pQ) is form-invariant under the action 
of SL(4,R), which constitutes the equivalence group of the problem. This action corresponds 
to linear transformations of the variables x,y,t and u. Since equivalence transformations pre- 
serve the integrability, any two SL(4, i?)-related equations are regarded as 'the same'. All our 
classification results are obtained modulo this equivalence. 

Based on the integrability conditions, in Sect. 4 we classify integrable equations of the form 
([T]) under various simplifying assumptions. This leads to a wide class of non-trivial examples, 
both known and new, which are expressible in terms of elementary functions, elliptic functions 
and modular forms. Just to mention a few of them, we have found an integrable equation 



UXX + Uyy ~ e Ut U tt = 0, 




a 



P'(U X ) ~ p'jlly) 
p(Ux)p(Uy) 



U x y + (3 



p'(u t ) - p'(u x ) 

p(u x )p(u t ) 



u x t + 7 



p'(Uy) - p'jUt) 
p{Uy)p{Ut) 



Uy t = 0, 
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here p is the Weierstrass p-function, (p') 2 = 4p 3 — g%, and a, (3, 7 are arbitrary constants. 
Another interesting example comes from the class 

u xy + (u x u y r(u t ))t = 0; 

for equations of this form the integrability conditions result in a single third order ODE for r, 

tit / 1 2\ 111 1 a 3 // 1 o /3 c 2 /2 n 

r [r — r ) — r + 4r r + 2r — or r =0, 

which appeared recently in the context of modular forms of level two pQ. Its generic solution is 
given by the Eisenstein series 

n=l 

which is associated with the congruence subgroup To (2) of the modular group. Further exam- 
ples of integrable equations expressible in terms of modular forms can be found in [22] in the 
classification of integrable Lagrangian equations of the form (pQ) corresponding to first order 
Lagrangian densities g{u u ,u y ,Ut). A generic Lagrangian density turns out to be a modular form 
of its arguments. Further generalizations of the above example include the equation 

u xy + 2(u x (log0)') t = 0, 

where 6 — ff) is the Jacobi theta-function, 

00 

6{z, t) = 1 + 2 e™ 2r cos(2vmz), 

71=1 

and prime denotes differentiation by u y . 

In Sect. 5 we study first order conservation laws, that is, relations of the form 

gi(u x ,Uy,U t ) x + g2(u X ,Uy,Ut)y + 9^{u x , % , U t )t = 0, (2) 

which hold identically modulo ([T]). Our second result (Theorem 2 of sect. 5) states that 

• Any integrable equation of the form ([1]) possesses exactly four first order con- 
servation laws. 

In Sect. 6 we investigate the existence of dispersionless Lax pairs (scalar pseudo-potentials), 

S t = F(S X , u x , u y , u t ), S y = G(S X , u x , u y , u t ), (3) 

which imply Eq. ([T]) via the consistency condition Sty = S y t (we point out that the dependence of 
F and G on S x is generally non- linear). Lax pairs of this type first appeared in the construction 
of the universal Whitham hierarchy, see [32^ [33] and references therein. It was observed in [H] 
that consistent Hamilton-Jacobi-type relations of the form © arise from the usual 'solitonic' 
Lax pairs in the dispersionless limit. Dispersionless Lax pairs constitute a key ingredient of the 
dispersionless 9-method, and a novel version of the inverse scattering transform [30l I35| . It 
was demonstrated in [18^120] that, for a number of particularly interesting classes of systems, the 
existence of dispersionless Lax pairs is equivalent to the existence of hydrodynamic reductions 
and, thus, to the integrability. Our third main result (Theorem 3 of Sect. 6) can be formulated 
as follows: 
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• Any integrable equation of the form (pQ) possesses a dispersionless Lax pair. 
Furthermore, the existence of a dispersionless Lax pair is equivalent to the 
existence of an infinity of hydrodynamic reductions and, thus, is necessary 
and sufficient for the integrability. 

Differential-geometric aspects of integrable equations of the form ([T]) are discussed in Sect. 7. 
Our main observation is that the equivalence group SL(4, R) acts by projective transformations 
on the space of first order derivatives p 1 = u x , p 2 = u y , p 3 = Ut, which is thus identified 
with the projective space P 3 , and the coefficient matrix fij(p) supplies P 3 with a well-defined 
conformal structure fijdp l dp 3 . The classification of integrable equations of the form ([1]) up 
to the action of the equivalence group SL(4,R) is therefore equivalent to the classification of 
conformal structures in P 3 up to the projective equivalence. Our first result in this connection 
is a characterization of linearizable equations (Theorem 4 of Sect. 7.1): 

• Equation of the form ([T]) is linearizable by a transformation from the equiva- 
lence group if and only if the corresponding conformal structure possesses a 
quadratic complex of null lines with the Segre symbol [(222)]. 

A tensorial characterization of linearizable and Lagrangian equations is provided in Sect. 7.2. 
The integrability conditions impose strong restrictions on the conformal structure fijdp l dp 3 
(Theorem 5 of Sect. 7.4): 

• Conformal structures associated with integrable equations of the form ([1]) are 
conformally flat. 

Thus, the theory of integrable equations of the form (pQ) has two 'flat' counterparts: one is the 
standard projective space P 3 with coordinates p 1 ,^ 2 ,^ 3 , and the projective action of SL(A,R); 
alternatively, one can speak of a projectively flat connection. Another is a flat conformal struc- 
ture fijdp l dpK Although, viewed separately, both structures are trivial, this is no longer true 
when they are imposed simultaneously: their 'flat coordinate systems' may not coincide (in 
fact, they do coincide for linearizable equations only). Our final result provides a differential- 
geometric characterization of the integrability conditions (Sect. 7.4): 

• Integrable equations of the form (pQ) correspond to conformal structures in P 3 
which possess an infinity of three-conjugate null coordinate systems parametrized 
by three arbitrary functions of one variable. 

Notice that we are in the realm of two different geometries, namely, conformal geometry (re- 
sponsible for the property of being null), and projective geometry (responsible for the property 
of being conjugate). It is quite remarkable that the moduli space of such structures is only 
20-dimensional! 

2 Derivation of the integrability conditions: the method of hy- 
drodynamic reductions 

As proposed in [T7] , the method of hydrodynamic reductions applies to quasilinear equations of 
the following general form: 

A(u)u x + B{u)u y + C(u)ut = 0; (4) 

here u = (it , u" 1 ) 1 is an m-component column vector of the dependent variables, and A, B, C 
are I x m matrices where I, the number of equations, is allowed to exceed the number of the 
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unknowns, m. The method of hydrodynamic reductions consists of seeking multi-phase solutions 
in the form 

u = u(R\...,R N ) (5) 

where the 'phases' R % {x,y,t) are required to satisfy a pair of consistent equations of hydrody- 
namic type, 

4 = fi l (R)K, R\ = X(R)K- (6) 

We recall that the consistency (or commutativity) conditions, R yt = R\ y , imply the following 
restrictions for the characteristic speeds fi l and X 1 : 

''hi' 1 _ 'V (7) 



y? - fi l \i - X 1 
i ¥= h di = d/drf, see 02]. 

Definition. |17j A system is said to be integrable if, for any number of phases N , 
it possesses infinitely many N -phase solutions parametrized by 2N arbitrary functions of one 
variable. 

Integrable equations of the form (j3J) arise in a whole range of physical and differential- 
geometric applications, and contain many particularly important classes of systems. Thus, in 
the case of square matrices, m = I, one arrives at systems of hydrodynamic type, 

u t + A{u)u x + B{u)u y = 0; 

we recall that n-phase solutions for hydrodynamic type systems, also known as solutions with a 
degenerate hodograph (when n = 1 or n = 2), have been extensively investigated in gas dynamics 
|43| . Second order quasilinear PDEs of the form (fT]) can be cast into the form @ by setting 
u = (u x ,u y ,ut); this gives a system of four equations in the three unknowns: m = 3, / = 4. 
Another interesting subclass corresponds to equations of the dispersionless Hirota type, 

FyUxx, U X y, Uyy, U X tl Uyt, Ufa) 0, 

see [21] and references therein: one can choose any five second order partial derivatives of the 
function u(x, y, t) as the dependent variables u. This corresponds to the case m = 5, / = 8. Let 
us illustrate the method of hydrodynamic reductions using the dKP equation as an example. 

Example. Rewriting the dKP equation, 

t^xt ^x^xx ^yy — 0, 

in the new variables a = u x , b = u y , c = ut, one obtains a quasilinear system of four equations 
in the three unknowns, 

a y = b x , a t = c x , b t = c y , a t - aa x - b y = 0. (8) 

Looking for iV-phase solutions in the form a = a(R 1 , R N ), b = b(R 1 , R N ), c = c(i? 1 , R N ) : 
where the phases R 1 satisfy Eqs. ©, one readily obtains the relations 

d i b = t M i d i a, d i c = X i d i a, A i = a + (/i i ) 2 . (9) 

The compatibility conditions didjb = djdib and didjC = djdiC imply 

dida = Ji^—d ia + S^dja, (io) 
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while the commutativity conditions (|7j) result in 

Ultimately, the substitution of (jlip into (jlOp implies the following system for a(R) and fJ- l (R), 



i 7^ J, which was first derived in |25|. I26j in the context of hydrodynamic reductions of Benney's 
moment equations. For any solution fx 1 , a of the system (I12p one can reconstruct A* and b, c 
by virtue of Q. A remarkable feature of the system (|12p is its multi-dimensional consistency: 
dkdj^i 1 = djdkfJ- 1 and dkdidja = djdidka. The general solution of the system (fl2l) depends, 
modulo reparametrizations R l — * f l {R l ) 1 on AT arbitrary functions of a single variable. Taking 
into account that N extra arbitrary functions come from the general solution of Eqs. ([6]), which 
can be obtained by the generalised hodograph method [44) . one ends up with an infinity of 
iV-phase solutions depending on 2N arbitrary functions of a single variable. We point out that 
the compatibility conditions, dkdjfi 1 = djdkfJ- 1 and dkdidjd = djdidka, involve triples of indices 
i/j/fc only. Thus, the consistency of the system (|12p for N = 3 implies its consistency for 
arbitrary N. This turns out to be a general phenomenon. 

The main result of this section is the following 
Theorem 1 The moduli space of integrable equations of the form HP is 20 -dimensional. 

Proof: 

Here we only sketch the proof: full details are provided in the Appendix. Our strategy is to 
derive a set of constraints which are necessary and sufficient for the existence of an infinity of 
re-phase solutions. These constraints constitute a complicated system of second order PDEs 
for the coefficients fy, which is in involution; after that, the calculation of the dimension of 
the moduli space reduces to a simple parameter count. The main steps of the derivation of 
the integrability conditions can be summarized as follows. First we introducing the variables 
a = u x , b = u y , c = ut, which transform Eq. (JT]) into the quasilinear form @, 

a y = b x , a t = c x , bt = Cy, fua x + f 2 2b y + /33c* + 2/i 2 o y + 2/130* + 2/ 23 6t = 0. 

Following the method of hydrodynamic reductions, we seek multi-phase solutions in the form 

a = a{R 1 ,...,R N ), b = b(R l , . . . , R N ), c = c(i2 1 , . . . , R N ), 

where the phases R 1 (x,y,t), R N (x,y,t) are arbitrary solutions of a pair of commuting 
hydrodynamic type flows (jHJ). The substitution of this ansatz into the above quasilinear system 
leads to the equations 

dib = [i l dia, die = \ l dia, 

along with the dispersion relation 

D(X i ,fi i ) = hi + /22O/) 2 + /asCA*) 2 + 2/ 12 // + 2/iaA* + 2/23//^ = 0. 
The consistency of the equations for dib and <9«c implies 
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Differentiating the dispersion relation with respect to R? , j ^ i, and keeping in mind Eqs. ([7]), 
one obtains explicit expressions for djX 1 and <9j/U* in the form 

djA* = (V - X' llij j;),,,. djfJ = (// - /r 1 ) !!,,(),<,. 

where B^ are certain rational expressions in A*, X J , fj, 1 , y? , whose coefficients depend on fij(a, b, c) 
and first order derivatives thereof (see the Appendix for explicit formulae). Thus, 

didja = —(Bij + Bji)di,adja. 

Calculating the consistency conditions djO^X 1 = dkdjX 1 , djO^fJ- 1 = dkdjfi 1 and didjd^a = 
didkdja, one arrives at 30 differential relations for the coefficients fij, which are linear in the sec- 
ond order derivatives thereof. These relations are manifestly conformally invariant, and without 
any loss of generality one can set, say, ^22 = 1- Solving for the second order derivatives of the 
remaining coefficients /n, /12, /13, /23, /33, one obtains 30 relations which can be represented in 
the symbolic form 

d 2 f ij = ^R(f kh df k i); (13) 

here F = det{fij}, and R is a polynomial expression quadratic in each of its arguments. In 
other words, any second order derivative of any coefficient fij is a certain explicit expression in 
terms of fki and first order derivatives thereof. A direct computation shows that the system 
(|13f) is in involution: all compatibility conditions are satisfied identically. Since the values of 
the five functions /n, /12, /13, /23, /33, and first order derivatives thereof, are not restricted by 
any additional constraints, we obtain a 5 + 3 • 5 = 20-dimensional moduli space of integrable 
equations. This finishes the proof of Theorem 1. 

The invariant tensorial formulation of the integrability conditions (|13l) is provided in Sect. 

7.4. 

Remark. Notice that, in two dimensions, any equation of the form 

fll(U x ,Uy)u xx + 2f 12 (u x ,Uy)u xy + h%{U x1 Uy)U y y = 

is automatically integrable. Indeed, in the new variables a = u x , b = u y it takes the form 
of a two-component quasilinear system, a y = b x , fn(a,b)a x + 2fi2(a,b)a y + f22(a,b)b y = 0, 
which linearises under the hodograph transformation interchanging dependent and independent 
variables. This trick, however, does not work in higher dimensions. 

3 The equivalence group 

The results of this section are not restricted to the dimension three, and hold in any dimension. 
Let us consider a multi-dimensional second order equation of the form 

/'/"<',<', = ( 14 ) 

where u = u(xi, ■■■,x n ) is a function of n independent variables, and the coefficients fij(p) = 
fij(p \ ■■■,p n ) depend on the first order derivatives p % = u Xi only. In matrix form, Eq. (|14p can 
be represented as 

trFU = 0, 

where F = (fij) is the (symmetric) matrix of coefficients, defined up to a scalar multiple, and 
U is the Hessian matrix of the function u. Our main observation is that the space P n with 
coordinates p 1 ,...,p n admits a natural projective action of the group SL(n + 1,R), and the 
matrix fij endows P n with a well-defined non-degenerate conformal structure fijdp l dp 1 . This 
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can be seen as follows. The class of equations (|14p is invariant under linear point transformations 
of the form 

x = Cx + bu, u = cx + flu, (15) 

where x = (x±, (column) vector of the independent variables, C is a re x n matrix, 
b and c are re-component column and row vectors, respectively, and (5 is a constant. The 
requirement that the transformation (I15p belongs to the special linear group SL{n + 1,R) is 
equivalent to the condition detC(/3 — cC b) = 1. The induced transformation law for the (row) 
vector of first order derivatives p = (p , ...,p n ) is manifestly projective, 

P = (pC-c)/(/3-p6), 
with the inverse transformation given by 

p = k pCT 1 + CC- 1 ; (16) 
here the scalar k is defined as k = f^pg-i^ = rf e t(c+bp) • ^ direct calculation gives 

# = Kdp (V 1 - ^p^J = *dp(C + ftp)" 1 . (17) 

Note also that dx = (C + 6p)dx. The transformation law for the Hessian matrix C7 can be 
obtained by substituting dp, dx into the equality dp = dx l U : 

U = -{C + bp) t U{C + bp). 

K 

Using this expression for U one obtains 

trFU = -tr[F(C + bpfUiC + ftp)] = -tr[(C + bp)F(C + bpW] = trFU 

K K 

where the new coefficient matrix F is given by 

F = -{C + bp)F(C + bp) t . (18) 

K 

Using the formulae (|17f) and ()18p one readily verifies the identity 

dpFdp 1 = KdpFdp 1 , 

which shows that the conformal class of the quadratic form dpFdp 1 = fijdp l dp 1 is defined in 
an invariant way. Thus, the coefficient matrix F(p) can be viewed as defining a conformal 
structure in the projective space P n with coordinates p 1 ,...,p n and the standard projective 
action of SL(n + 1, R). Thus, 

the study of equations of the form |j^[ ) up to linear transformations of the dependent and inde- 
pendent variables is equivalent to the study of conformal structures in P n up to the projective 
equivalence. The global counterpart of Eq. |j^[ ) would be a conformal structure on a projective 
manifold M n ( a manifold is said to be projective, or endowed with a flat projective structure, if 
it possesses an atlas with projective transition maps). 

The integrability conditions constitute a system of second order PDEs for the coefficient matrix 
F(p). Since point transformations preserve the integrability, the group generated by Eqs. (|16p 
and (|18p constitutes a point symmetry group of the integrability conditions. This SL{n + 1, R)- 
invariance plays a crucial role in the further analysis. In particular, all classification results 
presented below are formulated modulo this equivalence: two SL(n + 1, i?)-related equations are 
regarded as the 'same'. 

Returning to the case n = 3, we have a 20-dimensional moduli space of integrable equations, 
with the action of the equivalence group SL(4,R). One can prove that this action is locally 
free, that is, its generic orbits are 15-dimensional. Thus, up to the action of SL(4,R), a generic 
integrable equation is expected to depend on 5 essential parameters. 



8 



4 Examples and classification results 



Here we present some further examples (both known and new) and partial classification results 
of integrable PDEs of the form (JTJ) based on the integrability conditions derived in Sect. 2. 
Although these conditions are quite complicated in general, they can be written down and 
solved explicitly under various simplifying assumptions. A computer program which calculates 
the integrability conditions can be obtained from 

http : //www-staff . lboro . ac .uk/~maevf /bf t-supplementary-materials-2008 . tar .gz 

We construct a whole variety of new integrable examples expressible in elementary functions, 
elliptic functions and modular forms, thereby manifesting the richness of the problem. 



4.1 List of known examples 

In this subsection we bring together some of the known examples which arise, primarily, in the 
context of the dispersionless KP/Toda hierarchies. The simplest generalizations of the dKP 
equation are the modified dKP equation, 

1 2 

^xt ~i~ 2^x^j; ^yy UyU xx — 0, 
and the deformed dKP equation, 

i^xt ~\~ "^^x^xx tUyiixx ^yy u x u xx — 0, 

see e.g. [34J. Further examples include dispersionless limits of the KP and Toda singular 
manifold equations, 

Ut 3 / Uy \ \ 3 / Uy 



and 



u T I 4 V it 



Uxv — ~, Ti — ^Utt = U, 19) 

v (e n * - 1)(1 - e~ Ut ) v ; 

respectively [8j. Various symmetric forms of the dispersionless KP, BKP and Toda hierarchies 
were obtained in [HI [9] : 

(u y - u z )u yz + (u z - u x )u xz + (u x - u y )u xy = 0, 

f U , x(W"y ^z^jHyz Uy\V/% ^x^j^xz ^zi^x ^y^j^xy — O5 

(20) 

u x (u 2 - u 2 z )u yz + u y (u 2 z - u 2 x )u xz + u z (u 2 x - u 2 y )u xy = 0, 

(e n * - l)(e M * - e u *)u yz + (e u v - l)(e u * - e u *)u xz + (e Uz - l)(e u * - e u y)u xy = 0. 

A modification of the standard R-matrix scheme leads to the so-called r-th dispersionless mod- 
ified KP and r-Dym equations, 



3 — r (3 — r)(l — r) (3 — r)r (3 — r)(l — r) 2 

u xy ~ 7^ v? u tt ^ n utUxx + — ~ u x Uxt H u x u xx = 

(2 — r) z 2 — r 2 — r 



and 

3 - r / 1 1 



lit p. I n UyU XX U X U X y I 0, 

a 2 — r V 2 — r y 1 — r y ' 



9 



see [5], [6], [36]. Further examples arise in the context of the so-called 'universal hierarchy': 

- - ( }!±\ - ( 

Uyy — UyU X z U X yU z , U X z — U x Uyy U X yUy ; I I — I I j 

V U X J t \ U X J x 

see |38t 0U] . A class of integrable Euler-Lagrange equations of the form 

{g u Jx + {g Uy ) y + {g Ut )t = o 

was investigated in [20] ; here the Lagrangian density g is a function of the first order derivatives 
only, g = g(u x ,u y ,ut). Among the apparently new examples found in [20] one should mention 
the equation 

U t U X y + U X Uy t + UyU xt = 0, 

which corresponds to the Lagrangian density g = u x u y Ut- It was demonstrated in [22] that the 
generic integrable Lagrangian density g is a modular form of its arguments. 

Further examples of integrable equations of the form ([T]) result, via appropriate substitutions, 
from 2D integrable systems of hydrodynamic type. For instance, the paper [23] provides a 
complete list of integrable Hamiltonian systems of the form 

v t + (H v ) y = 0, w x + (H w ) y = (21) 

where the Hamiltonian density H is a function of v , w. Eq. (|2ip 9 implies the existence of a 
potential u such that Uy — w, u x — — H w . Expressing v and w in terms of u x and Uy and 
substituting these expressions into Eq. (|21|) i . one obtains an equation of the form ([T]). Thus, 
the simplest integrable example H = ^vw 2 results in the equation u x u y t — u y u xt + UyU yy = 0, 
etc. 

Another possible construction exploits the fact that any two-component integrable system 
of hydrodynamic type possesses a dispersionless Lax pair of the form 

S t = F(S x ,v, w), S y = G(S x ,v, w), 

see |18j . Here v and w are the dependent variables which satisfy a system of hydrodynamic type 
resulting from the compatibility conditions S y t = St y ■ Expressing v and w in terms of S x , S y 
and St, and substituting these expressions into the equations for v and w, one obtains a single 
second order equation of the form ([1]) for S. This construction is analogous to the 'eigenfunction 
equations' in soliton theory |29| . 

We have verified that all examples listed in this section indeed satisfy the integrability con- 
ditions of Sect. 2. 

4.2 Equations of the form u tt = f{u x , u y , Ut)(u xx + u yy ) 

For these equations, which can be viewed as nonlinear analogues of the (2 + l)-dimensional wave 
equation, the integrability conditions take the form 



Vfbb + fl ~ 




= o, 


2/ faa ~ 3/ a 


+ fS 


= o, 


2ff cc - 2// c 2 + f 2 a 


+f 2 b 


= o, 


fafc ~ 


facf 


= o, 


Zfafb - 


fabf 


= o, 


fbfc ~ 


fbcf 


= o, 
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recall that a = u x , b = u y , c = Uf. These relations are straightforward to solve. The case 
fa = h = leads, up to equivalence transformations, to a unique solution / = e c , which 
corresponds to the Boyer-Finley equation 



u t t = e Ut (u xx + u. 



yyi- 



Assuming f a and fb to be nonzero (notice that f a and fb can only be zero or nonzero simulta- 
neously), the last three equations imply 

f(a,b,c) - 



and the substitution into the remaining equations gives 

ae pc + ~/e~P c + 5 



f(a,b, c) 



a 2 + b 2 + £ 



where the constants a, f3, 7, 6, £ satisfy a single quadratic constraint 4q7 — 5 2 = 0. This corre- 
sponds to equations of the form 



ae Pu t + Je -f3u t + $ ^ 
ul + u 2 +i 



which are analogous to the singular manifold dToda equation (|19l) . 
4.3 Equations of the form u a = f(u x ,u y ,u t ) x 

These equations are related to the Hirota-type equations Qtt = Qxy, ^xt) via a substitu- 

tion u = £l x . In this form they were investigated in |41|, [T9] . The integrability conditions take 
the form 

n _ 2f cc , _ 2/ ac / cc ^fbcfcc 
Jccc — 7 ! Jacc — 7 > Jbcc 



fc fc fc 

f _ 2f ac r _ 2facfbc „ _ ^fbc 
Jaac — j, 1 Jabc — r i Jbbc — j, j 
Jc Jc Jc 

2 

fbbb = -J2 {fofbc + fbcifcfbb + 2/ac) ~ fcctfbfbb + 2/ a &)) , 
Jc 

2 2 

/af>6 = "72 (/a/be "I" fac(fcfbb ~l~ /ac) fcc{fafbb /aa)) j 
/c 

2 

/oab = "To \fcc\fbfaa ^fafab) fac(fbfac ^fcfab) fbcifcfaa 2/ a / ac )) , 
fc 

faaa = (if a fb )/ac fafbc fc (fab ~ faafbb) ~ fbbfeef a 

~^~facfc{faa ^(fafbb fbfab)) ~\~ 2fbc{fb(fcfaa fafac) fafefab) 
/cc((/o fb )faa ^fafbfab)^) ■ 
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This system is in involution, and its general solution depends on 10 arbitrary constants. The 
integration leads to the four essentially different canonical forms, 



' ■ — "2 , ri„-Au x 



f = u y + —(Au t + 2Bu x y + Ce 
^ u x ( u x Au x ) * u 2 1 Au x 



ut 6 



/ = In u y - In 6»i (u t , u x ) _ | / v( T )dr, 

see [41] . Here A, B, C are arbitrary constants, rj is a solution of the Chazy equation |12j . 

rj"' + 2rjrj" = 3(t/) 2 , 

and 9\ is the Jacobi theta-function. 



4.4 Equations of the form f(u x , u y , u t )u xy + g(u x , u y , u t )u xt + h(u x , u y , u t )u yt = 

Equations of this type arise in the context of the dispersionless KP, BKP and Toda hierarchies, 
see (|2(jp . Further examples were constructed in pj in the search for consistent triples of second 
order equations. First of all, the integr ability conditions imply that any equation within this 
class can be reduced to a simplified form, 

f(u x , u y )u xy + g(u x ,u t )u xt + h(u y ,u t )u yt = 0, 

via a multiplication by an appropriate scalar factor. In terms of the coefficients f(a,b),g(a,c) 
and h(b,c), the integrability conditions take the form 



Jab , J a t Jb i Ji 

h g g h 

_ h c f a f a h c 

9ac — ~~r9a H T~9c ~T~T~9i 

h f f h 

t 9c i . fb , fb 9c i 

he = —h + 7 /ic - -T—h, 
9 f f 9 

, _ gfa + f9a , ha ~ 9 fa f 9h - fhc , fgahc - 99ah , 
Jaa j J a + Jb ' Ja 1 ^ J 1 

, fhb + hf b fh b -hf b hg a - fg c , . fh b g c - hg a h b 

fob — 77 Jb H 7 Ja H 5 ^ 2u fo 

fh fg g 2 g 2 h 
gh c + hg c gh c - hg c hf a - gf b gf b h c - hf a h c 
g cc = — — g c H y — 9a H p 9c H Yh 9 > 

__ ha + gfa , gfa - ha , fh c -gh h , gf a h b -ff a h c 

9aa — y. 9a + ^ <?c + ^ 2 9a ~r 3) 

#/t c + %, - 0ft c 5/ 6 - /l/ a , , hf a g c ~ 9fb9c , 

n cc = r n c H 77 Ai6 H 75 n- c H 75 ra, 

gh fh f 2 f 2 g 

gh b + hf b hf b -fh b fg c - hg a hg a f b - ff b g c 

hbb = 77 h b H 7 h c H 5 ft & H —5 A. 

fh gh g 2 fg 2 

Remarkably, the first three equations appeared previously in [23] in the problem of classification 
of integrable Hamiltonian systems of hydrodynamic type. It was observed that their general 
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solution is representable in the form 



= p(a) - q(b) = r(c)-p{a) = q(b) - r(c) 
J P(a)Q(b) ' 9 P(a)R(c) ' Q(b)R(c) ' 1 ' 

where p, q, r and P, Q, R are arbitrary functions of the indicated arguments. The further analysis 
splits into four different cases depending on how many functions among p, q, r are constant. In 
the simplest case when all three of them are constant, the integrability conditions reduce to 
P" = Q" = R" = 0. Thus, any equation of the form 

a\u x u yt + a 2 u y u xt + a 3 u t u xy = 

is automatically integrable. Let us concentrate on the generic case when none of p, q, r are 
constant (intermediate cases can be considered in a similar way). Under the substitution (|22l) . 
the first three integrability conditions will be satisfied identically, while the last six imply a 
system of second order ODEs for p, q, r and P, Q, R: 



P" P'-Q' P'-R' Q'-R' 



p 


- q 


+ 


p 


— r 


q 


— r 


R' 


- Q' 


+ 


R' 


-P' 


Q' 


-P' 


r 


- q 


r 


~P 


q 


~P 


Q' 


-p' 


+ 


Q' 


-R' 


p' 


-R' 



(23) 



P^_ p _ q-P ( Rr' _ R ,\ + r-p / Qq' _ + p'P(q - 2p + r) 



pi g _ r y r _ p J r _ q yg _ p J ^_ r ^p_g) 
^ R = g-r / Pl_ _ A + p-r } Qqj _ A + r'R (p - 2r + q) 

R q — p \p — r ) p — q\q — r ) \ r — p)\ r — q) 

"Lq = p^l f Rr> _ #\ + r _zi (!£_ _ P '\ + q'Q(p ~ 2 i + r ) 



q' p — r \r — q J r — p \p — q J (q — p)(q — r) 

The separation of variables in Eqs. (|23p implies 

P' = F(p), Q' = F(q), R' = F(r), 

where F(-) is an arbitrary quadratic polynomial, F(x) = ex 2 + mix + mo; here the parameters 
e,m.i,mo play the role of separation constants. A similar separation of variables in Eqs. (|24|) 
results in 

p'P = S(p), q'Q = S(q), r'R = S(r), 

where £(•) is a cubic polynomial, S(x) = ex 3 + 712X 2 + n\x + uq. These equations lead to 
P'/P = p'F(p)/S(p), Q'/Q = q'F(q)/S(q), R'/R = r'F(r)/S(r), and the integration yields 

P = ^(p-aYip-Pfip—ff, Q = C2 (q- a r(q-f3y(q-^, R = c 3 (r - a^r - f3) u (X ~ if; 

here c, are three extra integration constants, a, f3, 7 are the roots of the polynomial S, and the 
exponents /x, v, 77, which are related to the coefficients mi, mo of the quadratic polynomial F, 
satisfy a single relation /i + v + r/ = 1. Thus, setting e = 1, we have 

p> = ±(p _ a )i-Ai(p _ p)i-»(p _ 7 )l-n, q' = l(q- a) l ~^(q - 0Y~ v {q - j) 1 ^, 
ci c 2 
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r' = -(r - af-^r - /3) 1 "^(r - 7) 1 -"; 
C3 

recall that quadratures of this type arise in the context of the Schwarz-Christoffel mappings 
of triangular domains. Particularly interesting examples correspond to the symmetric choice 
fj, = u = rj = 1/3. In this case the ODEs for p, q, r take the form 

P' 3 = \SHP), q' 3 = \s\q), r« = ±S 2 (r). ( 25 ) 
c l c 2 c 3 

We point out that the ansatz ([22]) possesses the obvious SL (2, R) -symmetry, 

ap + P aq + P ar + (3 _ P _ Q _ i? 
; F' <? ^ ; F> r ^ ; f> ^ ^ r- f> Q -> — f> ^ 



7P + 5 ' 75 + 5 ' 7r + (5 ' 71? + <5 ' 75 + 5 ' 7r + 5 ' 

which can be used to bring the polynomial S to a canonical form. In the case of three distinct 
roots one can reduce S to a quadratic, S(x) = x 2 + 53, so that the ODEs ([25]) imply p = 
p'(aiwi), (/ = p'{a2U2)i r = p'(asus) where 27af = 2/c 3 , and p is the Weierstrass p-function: 
(p') 2 = 4p 3 — 53 (notice that we are dealing with an equianharmonic case: gi = 0). This leads 
to the integrable equation 

p'(aiu x ) - p'{a2U y ) ^ | p'{azu t ) - p'{a\u x ) ^ | p^Qg%) - P^Qg^t) ^ _ Q 
p(aiw :c )p(a2'U2 / ) xy p(aiu :r )p(a3iii) x * p(a2U y )p(a 3 u t ) yt 

Up to an appropriate rescaling, this equation is equivalent to 



p'(u x ) - p'(uy) ^ | p P'jut) - p'(u x ) ^ ^ | p'(-^) - p'iut) 
p(u x )p(u y ) xy p(u x )p(u t ) xt p(u y )p(u t ) 



It possesses a degeneration g% — ► 0, p(x) — > 1/x 2 , resulting in 

au t (ul - Uy)u xy + (3u y (u 3 - u 3 x )u xt + ^u x (u 3 y - ul)u yt = 0, 
compare with ([20]) . 



4.5 Integrable equations in terms of modular forms and theta functions 

This section contains a number of more 'exotic' integrable examples which are not expressible 
in elementary functions. 

Let us begin with equations of the form 

(u y p{u t )) x + {u x q{u t )) y + (u x u y r(u t )) t = 0. 

The integrability conditions yield a complicated system of three third order ODEs for the func- 
tions p, q and r. Let us analyze special cases. 

Case 1. The choice q = p, r = p' corresponds to Lagrangian equations with the Lagrangian 
density u x u y p{ut). In this case the integrability conditions reduce to a single fourth order ODE 
for p, 

1111 1 2 11 r> /2\ 2 ///2 1 r> 1 11 III 1 o '3 /// n /2 //2 n 
p [p p — Ipp ) — p p + Ipp p p + op p — 9p p = U. 

It was shown in [22] that the generic solution of this equation is a modular form of level three, 
known as the Eisenstein series E\^{z): 

p(ut) = e (a2 ~ al3+l32)ut = 1 + 6e ut + 6e 3ut + 6e 4 "' + 12e 7 "' + 
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(the Eisenstein series E\^{z) results upon the substitution Ut — * 2iriz). It can also be written 
in the form 

/ e (3fc-l)« t e (3fc-2)w t \ 

p(u t ) = 1-6 I — e ( 3fc __ 1)Mt - 1 _ e ( 3k -2)ut J - 

Notice that a similar choice q = p,r = —p corresponds to equations of the form u xy —^-u x u y uu = 
0. Here the integrability conditions result in a fourth order ODE 

p 2 p"p"" - p 2 p"' 2 - 2pp" 3 + p l2 p" 2 = 0, 

whose properties are quite different from those of the above equation: first of all, one can reduce 
the order by setting p" = 2ph. This results in the second order ODE hh" — h' 2 — 2/i 3 = 0, which 
implies 

p 2su t 

h{u t ) = As 2 



s = const. This is the case of the singular manifold dToda equation (|19p . see also Ex. 2 of Sect. 
5. 

Case 2. Another interesting choice is p = 1, q = which corresponds to equations of the form 
u xy + {u x u y r(u t ))t = 0. The integrability conditions result in a single third order ODE for r, 

r"'(r' - r 2 ) - r" 2 + 4rV + 2r /3 - 6rV 2 = 0, 

which appeared recently in a different context in the theory of modular forms of level two: set 
r = y/2 to obtain Eq. (4.7) from [1]. This equation possesses a remarkable SL(2, i?)-invariance, 

az + (3 _ „ 2 , rN . _ 

z = -, r = Hz + o) r + 7(72: + 0), ad — 7/5 = 1; 

72: + 

here z = Ut- Modulo this SL(2, i?)-action, the generic solution is given by the series 

r(ut) = 1 + 8e 4ut - 8e Sut + 32e 12ut - 40e 16ut + 48e 20nt - 32e 24ut + 
which, upon setting e iut = q, coincides with the Eisenstein series, 

(-l) n nq n 



00 



n=l 

associated with the congruence subgroup To (2) of the modular group pp. 

Case 3. As a generalization of Case 2, let us consider equations of the form u xy +(u x r(u y , ut))t 
0. The integrability conditions take the form 

nbb _ rr bb - r bc + r 2 
r bb rr b - r c 

r bbc _ rr bc - r cc + r b r c 
r bb rr b - r c 

r b cc = ( 2rr bc ~ rr bbr cc ~ r cc r bc + r\r cc - 2r b r c r bc + 2r 2 r bb ), 

rr b - r c 

r C cc = ( 2r2r bc ~ 2r 2 r bb r cc + rr cc r bc + 4rrf r cc - 8rr b r c r bc + 

rr b -r c 

4rr 2 r bb - r 2 - r b r c r cc + 2r 2 r bc ). 
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Here the first two equations imply = a(rr^ — r c ), which is the well-known Burgers equation. 
Without any loss of generality we will assume a = 1. Under the substitution r = 2vb/v, 
the Burgers equation linearizes to the heat equation, v c = v^. Modulo this equation, the 
integrability conditions reduce to a single sixth order ODE for the function v, 

uVV 6 ) - uV 4 V 5 ) + 3v 3 v" 2 v^ - 5v 3 v"v"'v^ + 2v 3 v"' 3 - 3vWv^ - 2v 3 v'v"'v^ 
+ 5v 3 v'v^ 2 + 2v 2 v'v"v"' 2 + 6v 2 v' 2 v"v^ _ i 0u V VV 4 > + 2v W _ Gvv ' 2 v " 2 v "> 
+ 12vv' 3 v"' 2 - 6WV 5 > + 6v ,3 v" 3 - 12v'\"v'" + 6t/ V 4 > = 0, 

here prime denotes differentiation with respect to b. One can show that the generic solution 
of the heat equation constrained by this sixth order ODE is given by the formula v(u y ,ut) = 
9 (7?-, —^7) , where 6 is the Jacobi theta-function: 

00 

2, 



00, r) = 1 + 2 ^ e™ r cos(27rnz). 



71=1 



Remark. Further generalization, u xy + f(u x ,u y ,ut)t = 0, was discussed in [TH] in the disper- 
sionless Hirota form Q xy + f(Q x t, £lyt, &tt) = (these two representations are related via u = fit). 
It was demonstrated that the generic solution is given by a ratio of two Jacobi theta functions: 



f(u x ,u y ,u t ) = --In 



1 7 _ 9i(u t ,u y - u 2 



4 0l(Ut,Uy+U X ) 



5 Conservation laws 

Let us begin with some general remarks which clarify the geometric meaning of conservation 
laws of Eq. (|14p . Consider a first order conservation law in the form 

d Xl gi + ... + d Xn g n = 0, (26) 

where gi (p) are functions of the first order derivatives p 1 = u Xl , . . . , p n = u Xn . The requirement 
that the equality (f26|) is satisfied identically modulo Eq. ([Til) is equivalent to the set of relations 
d p igi = kfa, d p igj + d pJ gi = 2kfij, here k is a proportionality factor. These relations can be 
rewritten in the form 

dgxdp 1 + ... + dg n dp n = kf ij dp i dp>, (27) 

which has a clear geometric interpretation. Let us introduce a fiber bundle B with the base P n 
supplied with a flat conformal structure dgidp 1 (here p l are coordinates on the base P n , gi are 
coordinates along the fiber). The identity (|27p provides an embedding of the conformal structure 
fijdp i dp 3 , defined on the base, into the fiber bundle B with the flat conformal structure dgidp 1 . 
To uncover the geometry of B let us return to the equivalence group SL(n + 1,1?) which acts 
via (|15p on the space of independent variables X{. The induced action on conservation laws 
g = (01>->5n)* is 

g = k(C + 6p)g, 

recall that k = ^(g+fcp) ? see Sect. 3. Combining this with the transformation law 

dp = «;dp(C + bp)^ 1 

we obtain 

dp g = At 2 dp g. (28) 
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Introducing K = dp 1 A ... A dp n we can see that K = k u+1 K so that (|28p takes the form 

— 2 2 

dp g = dp g A"~^+i. (29) 

This shows that g has a natural interpretation as a 1-form with values in the canonical bundle 

2 „ . 2 

X "+ 1 , namely, g G H (P n , Q, <S) K Notice that for n = 3 (which is our main case of 

interest) this bundle appeared in [27] in the geometric theory of solutions of Einstein's equations. 

Let us return to the 3-dimensional case. Our main result is that any integrable equation of 
the form ([1]) possesses exactly four first order conservation laws of the form ([2]), 

gi{u x ,u y ,u t )x + g2{ux,Uy,u t ) y + gz{u x ,u y ,u t )t = 0. 

According to the discussion above, this means that there exists a 4-parameter linear family of 
sections g E ff°(P 3 , ft 1 (g) K~?) providing conformal embeddings of the conformal structure 
fijdp t dpi . Let us begin with illustrating examples. 

Example 1. The equation 

a\u x u yt + a 2 u y u xt + a 3 u t u xy = 0, 
which was found in Sect. 4.4, possesses four conservation laws of the form 

(o 2 + «3 - ai){u y u t ) x + (oi + a 3 - a 2 ){u x u t )y + (ai + a 2 - a 3 )(u x u y ) t = 0, 



a 2 [u t a 3 ii y + ai u t a 3 u x ) = 0, 



a l+ a 3 \ / a l+ a 3 

a 3 I u y a 2 u t ) + a x I My "2 n a . I = 0, 



a 2+ a 'i \ I a 2 + a 3 

a 3 I u x a i u t \ +a 2 [ u x a i -u, 



Example 2. The dToda singular manifold equation [5], 



-ntt = 0, 



^ (e n * - 1)(1 - e-"' 
possesses two conservation laws of the form 

as well as 

(p(Ut)Uy) x + {p(Ut)U X ) y - [p' \ut)U X Uy) t = 

where p satisfies the second order linear ODE 

p" e Ut 



2p (e u * - l) 2 ' 

as pointed out by M. Pavlov, it possesses two linearly independent solutions, 

e Ut + 1 , e Ut + 1 

p = and p = m — 2, 

y e"* - 1 e Ut - 1 
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which provide two extra conservation laws. 
Example 3. The equation 

ut(ul - v? y )u xy + u y (ul - u 3 x )u xt + u x (u y - Ut)Uyt = 0, 

which appeared in Sect. 4.4, possesses four conservation laws, 



(U X Uy(ul ~ Uy)) t + {U X U t {l4 - Ul))y + (UyU t (U y ~ Uf)) 3 



as well as 

Ut 



G 



u 



G 



U.r 



G 



u t 



G 



ut 



G 



Ut 



0. 



here the function G(s) is defined as G' = 3 ^_ 1 . Explicitly, one has 



here e = e 27ri//3 , 



G{s) = - (ln(s - 1) + e ln(s - e) + e 2 ln(s - e 2 )) 
1; notice that G is real- valued. 



Example 4. The equation 

p'{u x ) - p'(u y ) 



p (u t ) - p (u x ) p (u y ) - p Q t ) n 

p(u x )p[u^ * ' ■' 1 *• 1 ■' " ' 



"yj P(u x )p(u t ) "'" ' p(u 2/ )p(%) 

which appeared in Sect. 4.4, possesses four conservation laws, the first one being 

H(u y , Ut) x + H(ut, u x ) y + H(u x , Uy) t = 0, 
where the function H (r, s) is defined by the equations 

*2r.-\ , a, /„>C(*)p(a) - CMpM „ . 2f _. , , , , , , ,C(*)pO») - C(r)p( r ) 



flr = -rW-p(a)-4p(s)- 



p'(s) - p'(r) 



# s = C» + p(r)+4p(r)- 



p'(s) - p'(r) 



here £ is the Weierstrass zeta- function: £' = — p (we point out that the equations for H are 
automatically consistent: ff rjs = H sr ). Three extra conservation laws are of the form 

F(u x ,u y )t - F(u x ,u t ) y = 0, F(u y , U x ) t - F(u y , u t ) x = 0, F(u t ,u x ) y - F(u t , Uy) X = 0, 

where the function F(r, s) is defined by the equations 



Fr 



p(r)p(s) 



p'{s)-p'{r) 2 U ' 



p'(r) - p'(s)' 

Explicitly, one has (see [23], where the function F appeared in a different context): 

F(r, s) = - (in air — s) + e In aier — s) + e 2 In a(e 2 r — s)) ; 
6 

here e 3 = 1 and a is the Weierstrass sigma-function: a' /a = £. Notice that F is real- valued. 
The main result of this section is the following 

Theorem 2 Any integrable quasilinear PDE of the form ([T]) possesses four first order conser- 
vation laws. 
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Remark 1. We would like to stress that the converse statement is not true: the existence of 
four conservation laws does not necessarily imply the integrability. For instance, any Euler- 
Lagrange equation of the form {g Ux )x + {du y )y + (9u t )t = 0, corresponding to the Lagrangian 
density g(u x ,u y ,ut), is manifestly conservative and, moreover, possesses three extra first order 
conservation laws 

{u x g Ux ~ g)x + (u x g Uy ) y + (u x g Ut ) t = 0, 
(u y g Ux )x + {u y g Uy - g) y + {u y g Ut )t = 0, 

(u t g Ux )x + (utg Uy ) y + (u t g Ut - g) t = 0, 

which constitute components of the energy-momentum tensor. On the other hand, as shown 
in [20| . the integrability conditions are very restrictive and reduce to a complicated system of 
fourth order PDEs for the Lagrangian density g, resulting in the finite-dimensionality of the 
moduli space of integrable Lagrangians. 

Remark 2. Integrable multi-dimensional equations normally possess infinite hierarchies of 
higher order conservation laws, which are generically non-local. The discussion of higher con- 
servation laws is beyond the scope of our paper. 

Proof of Theorem 2: 

Differentiating Eq. ([2]) we obtain 

(gi) a u xx + {g2)bu yy + (gz) c u zz + ((gi)b + {g2) a )u xy + ((gi) c + (g3) a )u xz + {{g 2 ) c + (53)5)%*: = 0. 
Let us introduce the new variables s l ,r l via 

r 1 = (gi) a , r 2 = (g 2 ) b , r 3 = (g 3 ) c , 

2s 1 = (g 2 ) c + (g 3 )b, 2s 2 = (gi) c + (g 3 ) a , 2s 3 = ( 5l ) 6 + {g 2 ) a ; 
one can verify that these variables automatically satisfy the compatibility conditions 

s aa = s ab + s ac ~ r bci ^ s ab = r bb + r aa' 

s bb = s ab + s bc ~ r "aci ^ s bc = r cc + r bb' (30) 
^cc $ ac + Sfo c r abi 2s ac r cc + T aa . 

For a relation of the from ([2]) to be a conservation law of Eq. (fTJ) one has to require the existence 
of a factor k(a,b, c) such that 

r 1 = kfn, r 2 = fc/22, r 3 = fc/ 33 , 
s 1 = kf 23, s 2 = kf 13 , s 3 = kf 12 . 

Substituting this into the compatibility conditions (|30|) we arrive at six linear equations for 
the 'integrating factor' k. Solving this system for the second order partial derivatives kij (the 
variables a,b,c are labelled by indices 1,2,3, respectively), one obtains 

kij = —(fisfrj + fjsfri ~ fijfrs)(kfpq,M + 2k p ftl t q)c P e 9 . (31) 

Here e l i k is the totally antisymmetric tensor dual to the volume form of the metric corre- 
sponding to the equation ([T]), that is, e 123 = l/y/F, e 213 = —I/VP, etc, F = det/y. The system 
(|3ip appears to be in involution modulo the integrability conditions (|13p . Since the variable k 
and its first order derivatives k a ,kb,k c are not restricted by any additional constraints, there is 
a 4-parameter freedom for the integrating factor. This finishes the proof. 
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Example 5. Let us write out Eqs. (|3ip for the dispersionless Hirota equation discussed in Ex. 
1. We have fa = fa = /33 = 0, f\ 2 = 03c, /13 = 026, fa = a\a, and the substitution into (f3T1) 
implies 

3 3 3 

kab = k ac = kb c = 0, k aa = k a , kbb = —ykfj, k cc = k c . 

a c 

This leads to the four values for the integrating factor: k = const, k = a -2 , k = b~ 2 , k = c -2 , 
which correspond to the four conservation laws from Ex. 1. 



6 Dispersionless Lax pairs 

In this section we prove that any integrable equation of the form (pQ) possesses a dispersionless 
Lax pair, 

S t = F(S X , u x , u y , u t ), S y = G(S X , u x , u y , u t ). 
Let us begin with illustrating examples. 
Example 1. The dispersionless Hirota equation, 

a\u x u yt + a 2 UyU xt + a 3 u t u xy = 0, a% + a 2 + a 3 = 0, 

possesses the Lax pair 

St S X Sy S X ^ 

U t U X Uy U X 

here the constants A and fx satisfy a single quadratic constraint ai\/i + a 2 [i + a 3 X = 0. Notice 
that this Lax pair is linear in S. It was discussed in |46j in the context of Veronese webs, and was 
used to solve the dispersionless Hirota equation via a non-linear Riemann problem. We will see 
below that the case a\ + a 2 + a 3 ^ is far more complicated, leading to Lax pairs parametrized 
by hypergeometric functions. 



Example 2. The dToda singular manifold equation, 



u r u 



x 



-u a = 0, 



y (e n * - 1)(1 - e-«*) 
possesses the Lax pair 

St = F(S x /u x ,u t ), S y = F 2 (S x /u x ,u t )u y ; 
here the function ^(x 1 ,^ 2 ) is defined as 

F(x\x 2 ) = -- In - Ax x (l - e" x2 )) , A = const, 

F 2 = dF/dx 2 . Although the constant A can be eliminated by a rescaling S — ► XS, one can 
consider the limit as A — ► 0. This results in a linear Lax pair for the same equation: 

5, = (l-e-^)-, Sy = e- U ^u y . 

Parametric Lax pairs. In many cases it turns out to be more convenient to work with 
parametric Lax pairs, 

S x = f(p,u x ,u y ,u t ), S y = g(p,u x ,u y ,ut), S t = h(p,u x ,u y ,ut), (32) 

here p is a parameter. Expressing p from the first equation and substituting into the last two one 
gets a Lax pair in the form ([3]). Similar parametric Lax pairs appeared previously in the context 



20 



of the universal Whitham hierarchy [32], see also |39| . The condition for (|32p to be a Lax pair 
for Eq. (JXJ) can be derived as follows. Thinking of p as a function of x, y, t, and calculating the 
compatibility conditions S xy = S yx , S x t = St x , S y t = Sty, one obtains three relations which are 
linear in p x ,p y and pt. It is easy to see that the rank of the matrix at the derivatives p x ,p y ,Pt 
equals two, so that one can obtain a single relation which does not contain the derivatives of p: 

(h a g P ~ g a h p )u xx + (f b h p - h b fp)u yy + (g c f p - f c gp)u u + 
((fa ~ 9b)h p + h b g p - h a f p )u xy + 

(33) 

((h c - f a )g P + g a f P - g c hp)u xt + 

((gb - h c )fp + f c h p - f b gp)u y t = 0; 
this relation must be satisfied identically modulo Eq. ([T]). This requirement leads to the relations 

Kg P - g a h p = kfn, f b hp- h b f p = kf 22 , g c f P - f c g P = fc/33, 

(fa ~ g )h p + h b g p - h a f p = 2fe/i2, 

(34) 

(h c - fa)g P + g a f P - g c h p = 2kf 13 , 

(gb - h c )fp + f c hp - f b g p = 2kf 2 3, 

where k = k(p,a,b,c) is the coefficient of proportionality. The set of relations (|34j) can be 
represented in a compact form 

/ f P g P hp \ 
det df dg dh = kfijdp L dp ) , 
\ da db dc J 

recall that (u x ,u yi u t ) = (a,b,c) = (p 1 ,p 2 ,p^)- We point out that, by virtue of ([341 . the triple 
(fp,g p ,h p ) satisfies the dispersion relation: 

fnfp + /225p + /33/ip + 2fi 2 f P g P + 2f 13 fph p + 2f 23 g p h p = 0. 

Parametric Lax pairs are particularly useful when the equation under study is symmetric under 
the interchange of x, y, t: 

Example 3. Let us consider the equation 

a\u x Uy t + a 2 u y u xt + a 3 u t u xy = 0, 

where, in contrast to Ex. 1, the constants a±, a 2 , 03 are arbitrary. Without any loss of generality 
we will normalize them so that a\ + a 2 + a 3 = 1. We seek a Lax pair in parametric form 

S x = f(p)u x , S y = g(p)u y , S t = h(p)u t . 

The corresponding Eq. ([33]) is 

(g - h)f'u x u yt + (h- f)g'u y u xt + (f - g)h'u t u xy = 0, 

here ' = d/dp, so that one can set 

f' = ai k(f-g)(f-h), g ' = a 2 k(g-f)(g-h), ti = a 3 k(h - f)(h - g), (35) 
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where k = k{p) is a coefficient of proportionality. We point out that the system f)35[) possesses 
a conservation law 

( g -h)^{h-f)^{f-g)^= const. 

Introducing the new independent variable q = (f—g) / {f—h) (notice that we have a reparametriza- 
tion freedom p — > (p(p)), and using the identity g' = &(/ — — h)/{f — h), one can linearize 
the system (|3"5j) . 

3-=i -i- = a 2(h- f), — = — (f-h), 

dq I — q dq dq q 

Noticing that f — h= (l — q)~ ai q~ a3 , one can reduce these equations to hypergeometric integrals, 

?L = a l {l-q)- a ^ 1 q- a \ f = -a 2 (l - q^q^, ^ = o 3 (l - g)- 01 ^- 1 . 
Explicitly, one arrives at the 'reparametrized' Lax pair 

- = /(«), ^ = <?(<?), - = %), (36) 

Mx My Mj 

where 

f = F(q), g = F{q)-{l-qy^q l - a \ h = F(q) - (1 - q)-^q-^ , 
and F(q) is a the hypergeometric function: 

F(0) = T^9 1_a3 2^1(01 + l,l-a 3 ; 2 -a 3 ;g). 
1 - a 3 

In the symmetric case a\ = 02 = a 3 = \ we obtain a Lagrangian equation 

+ %M x t + u t u xy = 0, 

which corresponds to the Lagrangian density u x u y ut- Its parametric Lax pair was calculated in 
|20j in the form (here the parameter is denoted by z), 

^- <W , 3l_ <W + £M+*, £_ (W + £!Mz* (37) 
Mx My 2p(z) Mi 2p(z) 

where p(z) is the Weierstrass p-function, (p 1 ) 2 = + A 2 (notice that g 2 = 0, 53 = —A 2 ), and 
((z) is the corresponding zeta- function: Q' = —p. The Lax pair (|37|) transforms into ([3"o]) via a 
substitution g = (p'(-z) + \)/{p'{z) — A). 

Example 4. The equation 

m 4 (m^ - Uy)u xy + m 2/ (m| - M^)u x . t + u x (u y - ul)u yt = 
possesses the parametric Lax pair 

s x = g(^), s v = g(<], s t = G ( u -i 



p j \ p j \p 

where the function G(s) satisfies the equations G' = -rzr- Explicitly, one has 

G(s) = - (ln(s - 1) + e ln(s - e) + e 2 ln(s - e 2 )) . 
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Example 5. The equation 



'{Ux) ~ p'(Uy) p'(Ut) ~ P'(U X ) p'(Uy) ~ p'(Ut) 



p{u x )p{Uy) ■ p{u x )p{Ut) ' p{Uy)p{Ut) 

possesses the parametric Lax pair 

S x = F(p,u x ), Sy = F(p,u y ), S t = F(p,u t ), 
where the function F(r, s) is defined by the equations 



= p(r)p(s) = P 2 (r) _ 1 

r p'(r)-p'( S y 3 p'(s)-p'(r) 2 U ] 



Explicitly, one has 



Fir, s) = - (in air — s) + e In aier — s) + e 2 In a(e 2 r — s)) , 
6 

compare with Ex. 4 of Sect. 5. Similar parametric Lax pairs can be constructed for all integrablc 
examples obtained in Sect. 4.4. 

The main result of this Section is the following 

Theorem 3 Any integrable equation of the form (pQ) possesses a dispersionless Lax pair. Fur- 
thermore, the existence of a dispersionless Lax pair is equivalent to the existence of an infinity 
of hydrodynamic reductions and, thus, is necessary and sufficient for the integrability. 

Proof: 

Setting S x = £, u x = a, u y = b, ut = c, and calculating the consistency condition for the Lax 
pair ([3]), one obtains 

Sty — S y t = (F^G a — G^F a )u xx + F b u yy — G c uu+ 

(F a + F^Gb - G^F b )u xy - {G a + Gt.F c - F^G c )u xt + (F c - G b )u yt . 
Since this expression has to vanish modulo Eq. (pQ), one arrives at the relations 

F i G a -G^F a = kf n , F b = kf 22 , G c = -kf 33 , 

F a + F^G b -G^F b = 2kf 12 , G a + G^F c -F^G c = -2kf 13 , F c -G b = 2kf 23 , 

where k = k(^, a, b, c) is the coefficient of proportionality. The last five relations imply the 
expressions for the derivatives of F and G in the form 

F a = 2kf 12 + G^kf 22 + Fs(kf 23 - p), G a = -2kf 13 - F^kf 33 - G(.{kf 23 +p), 

F b = kf 22 , Gb = -kf 23+p , (38) 

F c = kf 23 +p, G c = -kf 33 , 

where p = p(£, a, b, c) is yet another auxiliary function. Substituting these expressions into the 
first relation, one can see that F^ and G^ have to satisfy the dispersion relation, 

hi + h2G\ + f 33 Fl + 2/ 12 G e + 2f 13 F ( + 2f 23 F^ = 0. (39) 

To close the system (|3"5|) - (|3"9"|) one proceeds as follows. Calculating the consistency conditions 
for Eqs. (|38j) . F ab = F ba , G a b = G ba , etc, six conditions altogether, and differentiating the 
dispersion relation (|39p by a, b, c and £, one obtains ten relations which can be solved for F^, G^ 
and the first order derivatives of k and p. The resulting system is in involution if and only if the 
integrability conditions of Sect. 2 are satisfied. This finishes the proof of Theorem 3. 
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7 Differential-geometric aspects of the integrability conditions 



As explained in Sect. 3, the differential-geometric picture behind equations of the form (|14H is the 
projective space P n with coordinates p , ~.,p n supplied with the conformal structure fijdp l dpi . 
The equivalence group SL(n + 1,R) acts by projective transformations of P n . Two equations 
are equivalent if and only if conformal classes of the corresponding metrics are projectively 
equivalent. 

Let us consider Lagrangian equations of the form (I14p which arise as the Euler-Lagrange 
equations from the functionals f g(p)dx where the density g(p 1 , ...,p n ) depends on the first 
order derivatives p l = u Xi only. In this case the coefficient matrix /»,■ is the Hessian matrix of 
g. Our first remark, which is true in any dimension, is that the class of Lagrangian systems is 
invariant under the action of the equivalence group defined by Eqs. (|16p . (|18j) . One can show 
that the extension of the projective action (|16p to the Lagrangian density g is given by the 
formula 

' = TTiW (40) 

The projective invariance of the class of Lagrangian systems can be seen as follows. Let us 
consider the Lagrangian density g(p) as the equation of a hypersurface in P n+1 defined as 
p n+1 = g(p). The second fundamental form of this hypersurface coincides with the conformal 
class of the second differential d 2 g. The transformation (|16p . (|40p is a projective transformation 
in P n+1 . Thus, the fact that d 2 g transforms into d 2 g (up to a conformal factor) is nothing but 
the well-known projective invariance of the second fundamental form. 

A geometric characterization of linearizable equations (|14p is provided in Sect. 7.1 (Theo- 
rem 4). To be precise, we will be interested in those equations which can be linearized by a 
transformation from the equivalence group. 

A simple tensorial characterization of linearizable and Lagrangian equations is given in Sect. 
7.2. The answer is formulated in terms of the tensor and the projectively flat connection V 
with ChristoffePs symbols Tj k = Sj5 l k + Sfc<5*- which are naturally assocated with the conformal 
structure fijdp l dp J . 

Invariant differential-geometric formulation of the integrability conditions (|13|) derived in 
Sect. 2 is provided in Sect. 7.3. This involves the tensor a^j- and its covariant derivative with 
respect to V. 

Finally, a simple differential-geometric characterization of conformal structures correspond- 
ing to integrable equations is proposed in Sect. 7.4. 

7.1 Linearizable equations and quadratic line complexes 

Before formulating the main result, let us summarize the properties of linear (linearizable) 
equations. 

Example 1. Consider the 3-dimensional linear wave equation, 

notice first that it is Lagrangian with the quadratic Lagrangian density g = u 2 + u 2 — u 2 . The 
associated conformal structure in P 3 corresponds to the standard Lorentzian metric (dp 1 ) 2 + 
(dp 2 ) 2 — (dp 3 ) 2 ; here p 1 = u x , p 2 = u y , p 3 = ut- This conformal structure possesses a 3- 
parameter family of null lines defined by the equations 

p 1 = ap 3 + (3, p 2 = 7p 3 + 5, (41) 
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where the constants a, f3, 7, 6 satisfy a single quadratic constraint a 2 + 7 2 = 1. The last property 
can be reformulated in a projectively invariant way as follows. Recall that a 3-parameter family 
of lines in P 3 is called a line complex. A complex is said to be quadratic if it is defined by a 
single quadratic relation among the Plucker coordinates in the space of lines (which is identified 
with the Plucker quadric in P 5 ). In the parametrization (|4ip . the Plucker coordinates are 

(1 : a : (3 : 7 : 5 : aS — 7/3). 

Fixing a point in P 3 with coordinates pj, p 2 ,, Pq, the lines of the complex passing through this 
point generate a quadratic cone with the equation (p 1 — pj) 2 + (p 2 — p 2 ,) 2 = (p 3 — pjj) 2 ; these 
cones are nothing but the null cones of the corresponding conformal structure. Introducing in P 3 
homogeneous coordinates q° : q 1 : q 2 : q 3 viap* = q l /q°, one can see that the intersection of null 
cones with the plane at infinity (defined by the equation q° = 0) is the conic (q 1 ) 2 + (q 2 ) 2 = (<Z 3 ) 2 - 
Thus, all quadratic cones of our complex pass through one and the same plane conic (complexes 
of this type have the Segre symbol [(222)], see [28ll3j). Summarizing, we see that 

(a) the linear wave equation is Lagrangian; 

(b) the corresponding conformal structure possesses a three-parameter family of null lines which 
form a quadratic complex with the Segre symbol [(222)] (equivalently, quadratic cones of the 
complex pass through one and the same plane conic). 

Reformulated in these terms, both properties are manifestly projectively invariant, and hold for 
arbitrary equations related to the linear wave equation via the action of the equivalence group. 
Example 2. Let us consider the equation 

(1 - It 2 - Uy)u tt - U\(U XX + Uyy) + 2U t (U X U xt + UyUyt) = 0, 

u 2 +u 2 — 1 

which is Lagrangian with the Lagrangian density g = - u ^ — . The corresponding conformal 
structure, 

(1 - (p 1 ) 2 - {p 2 ) 2 )(dp 3 f - {p^fiidp 1 ) 2 + (dp 2 ) 2 ) + 2p 3 {p 1 dp 1 dp 3 + p 2 dp 2 dp 3 ), 

possesses a 3-parameter family of null lines (I4ip specified by a single quadratic relation /3 2 + 5 2 = 
1. This defines a quadratic complex whose null cones pass through one and the same planar 
conic defined by the equation (q 1 ) 2 + (q 2 ) 2 = (q ) 2 in the plane q 3 = (recall that q % are 
homogeneous coordinates in P 3 ). The equation linearizes (to the wave equation from Ex. 1) 
under the transformation 

u = —t, t = —u, y = y, x = x, 
which generates a projective transformation of the derivatives, 

1 ~ ~ *a x 

Uf , Uy , U X . 

u t u t u t 

This extends to the transformation of the Lagrangian densities as g = g/ut- One can verify 
that the quadratic Lagrangian density g of the linear wave equation transforms to the density 
g of the linearizable equation from Ex. 2. Geometrically, this transformation is nothing but a 
projective transformation which sends the plane q 3 = to the plane at infinity. 

Our observations are summarized in the following theorem which, in fact, holds in any 
dimension. 

Theorem 4 The following conditions are equivalent: 

(1) Eq. is linearizable by a transformation from the equivalence group. 
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(2) Eq. |7^p is Lagrangian with the Lagrangian density g = j^y where Q and I are arbitrary 
quadratic and linear forms in p 1 , ...,p n , respectively (not necessarily homogeneous). 

(3) The conformal structure fijdp t dp J possesses a complex (that is, a (2n — 3) -parameter family) 
of null lines whose quadratic cones pass through a stationary hyperplane quadric. For n = 3 
these conditions are equiavalent to the requirement that the complex has Segre symbol [(222)]. 

Proof: 

The equivalence of (1) and (2) can be seen as follows. Suppose that Eq. (fT4"|) is linearizable. 
Then the corresponding conformal structure fijdp l dp 3 is transformable to a constant coefficient 
form. Since any constant coefficient linear system is Lagrangian with a quadratic Lagrangian 
density g, and the class of Lagrangian systems is projectively invariant, any linearizable equation 
is necessarily Lagrangian. Applying the projective transformation (|16p . (|40p to a quadratic 
Lagrangian density, one obtains a density of the form g = where Q and I are quadratic 
and linear expressions in p 1 ,...,p n , respectively. Conversely, given a Lagrangian density of the 
form g = ®Pw , and applying any projective transformation which has the linear form Z(p) in 
the denominator, one obtains a purely quadratic Lagrangian density which gives rise to a linear 
equation. This establishes the equivalence of (1) and (2). 

The implication (1) ==> (3) is straightforward: any constant coefficient conformal struc- 
ture possesses a complex of null lines whose quadratic cones pass through a stationary quadric 
belonging to the hyperplane at infinity. Conversely, consider an equation whose conformal struc- 
ture possesses a quadratic complex of null lines such that all null cones pass through a stationary 
quadric belonging to a stationary hyperplane H. Applying a projective transformation which 
sends H to a hyperplane at infinity, we obtain a linear equation with constant coefficients. In the 
case n = 3 one can also refer to the Proposition 4.3. of [I] which implies that, up to projective 
equivalence, there exists a unique quadratic complex with the Serge symbol [(222)]. 

The constraints on the complex are crucial for the linearizability. 
Example 3. Let us consider the dispersionless Hirota equation 

a\u x u yt + a 2 u y u xt + a 3 u t u xy = 0, a\ + a 2 + a 3 = 0; 

we point out that the equation is integrable for any values of constants, not necessarily satisfying 
this relation. The corresponding conformal structure 

a±p dp dp + a 2 p dp dp + a 3 p dp dp 

possesses a three-parameter family of null lines which form a quadratic complex defined by the 
equation a\(3^ + a 2 a5 = 0. This complex is not of the Segre type [(222)], therefore, the equation 
is not linearizable. In fact, it is easy to see that it is not Lagrangian. 
Example 4. Let us consider the equation 

UXX + (Uy - l)U U - 2UyU t Uy t + ufUyy = 0] 

one can show that it is not integrable, and not Lagrangian. The corresponding conformal 
structure 

{dp 1 ) 2 + ((p 2 ) 2 - l)(dp 3 ) 2 - 2p 2 p 3 dp 2 dp 3 + (p 3 ) 2 {dp 2 ) 2 

possesses a three-parameter family of null lines which form a quadratic complex defined by the 
equation a 2 + 5 2 = 1, which is again not of the Segre type [(222)]. 

Remark. The condition for a conformal structure fijdp l dp 3 in P n to possess a complex of null 
lines is equivalent to a simple differential-geometric constraint 

9 (kfij) = <P(kfij)l ( 42 ) 
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here d k = d p k, <p k is a covector, and brackets denote a complete symmetrization in i,j,k. 
Contracting (|42p with f lJ we obtain 

<Pk = ^^f pq (d k f pq + 2d p f gk ). 

The condition (|42l) characterizes conformal structures coming from quadratic line complexes, see 
e.g. O H2] and references therein. Thus, the identity (|42p is necessary (although not sufficient) 
for the linearizability. An invariant characterization of linearizable equations is provided in Sect. 
7.2. below. 



7.2 Tensorial characterization of linearizable and Lagrangian equations 

The results of this section are valid in any dimension, and provide a simple differential-geometric 
criterion of the linearizability for equations of the form (|14|) . 

Lemma. An equation of the form is linearizable by a transformation from the equivalence 
group SL(n + 1, R) iff there exists a flat connection V with Christoffel symbols Y l - k = Sj5 k + s k S l j 
such that 

^kfij = Ckfij] (43) 
(connections satisfying Eq. |J5|) are known as Weyl connections). 

Proof: 

The necessity is straightforward: given a linear equation with constant coefficients fij, we im- 
mediately arrive at (|43p where V is a flat connection with zero Christoffel symbols, = d k = 
d/dpk, and cj~ = 0. Applying a transformation from the equivalence group (which acts projec- 
tively on the space P n with coordinates p 1 , ...,p n ) to a flat connection V, we will obtain a flat 
connection with nonzero Christoffel symbols of the form T l - k = Sj5 k + Skfij, which still satisfies 
Eqs. (|43p . Moreover, the condition (|43p is manifestly invariant under rescalings fij — ► rfij 
(under such rescalings, the covector c k transforms to + Vfelnr). 

Conversely, suppose a connection V has Christoffel symbols of the form T l - k = Sj5 k + SkSj 
(connections of this form are known as projectively flat: their geodesies are straight lines). If, in 
addition, V is flat (has zero curvature tensor), there exists a projective transformation bringing 
Christoffel symbols to zero. In the new coordinates Eq. (|43p will take the form d^fij = c^fij, 
which implies that the coefficient matrix is proportional to a constant matrix. This finishes 
the proof. 

Relations (I43p lead to explicit tensorial constraints for as follows. Taking into account 
that T l - k = SjS k + Sk5j one can rewrite (@3|) as 

dkfij = (cfe + 2s k )fij + Sif kj + Sjf ki . (44) 

Contacting Eqs. (|44|) with the inverse matrix f lJ one arrives at the relations 

f 3 d k fij = nc k + 2(n + l)s k , f j djf ik = c k + {n + 3)s k , 

with a double summation over i and j. This implies 

ndjfik) , 

(45) 

2(n + l)djf lk ). 



Sk ~ (n+2)(l-n) (®kfij 
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Given an arbitrary conformal structure fijdp l dp J in P n , let us introduce the tensor 

a ijk = dkfij ipk ^kjfij 'Si/fcj Sjfkii 

here s^, are the same as in Eqs. P5j) . One can readily verify the apolarity relations f l ^&ijk — 
0, f lJ aikj = 0. Thus, we can formulate the following 

Proposition 1. Equation |1^[ ) is linearizable by a transformation from the equivalence group 
SL{n + 1,R) iff the corresponding conformal structure satisfies the following two properties: 

(1) the tensor a^k vanishes; 

(2) the connection T l - k = Sj5 l k + Sfc<5* is flat; a simple calculation shows that this condition is 
equivalent to djSi — SiSj = 0. 

Remark 1. In a somewhat different form, the tensor a^k appeared previously in [12] in the 
study of manifolds of quadratic cones in P n . It was proved that the vanishing of alone 
implies the existence of a stationary hyperquadric Q n_1 C P n such that all cones of the family 
are tangential to Q n . Clearly, all such conformal structures are projectively equivalent, and 
can be brought to a canonical form 

[(p,p)-l](dp,dp)-(p,(ip) 2 = 0, 

where p = (p , ...,p n ) are coordinates in P n , and ( , ) is the standard scalar product. The null 
cones of this conformal structure are tangential to a unit sphere centered at the origin. The 
corresponding second order equation takes the form 

[(Vu) 2 - l}Au - (Vu)F(Vu)* = 0, 

where Vu = (u x i, ...,u x n) is the gradient of u, A is the Laplacian, and H is the Hessian matrix 
of u. In the three-dimensional case we arrive at the equation 

(u\ + U 2 - 1)U XX + (U 2 X + U 2 - l)u yy + (U 2 X + U 2 y - \)utt 

-2(U X UyU X y + U X U t U xt + UyU t Uyt) = 0] 

notice that it is Lagrangian: the corresponding Lagrangian density yjl — u x — Uy — u 2 governs 

minimal hypersurfaces z = u(x,y,t) in the Lorentzian space with the metric dx 2 +dy 2 -\-dt 2 — dz 2 . 
We have verified that this equation does not satisfy the integrability conditions of Sect. 2. 

Remark 2. Another result of [12] (also formulated in different terms) states that, imposed si- 
multaneously, conditions (1) and (2) imply the existence of a stationary hyperplane H in P n with 
a stationary quadric Q n ~ 2 C H such that all cones of the family pass through Q n ~ 2 . This pro- 
vides an alternative projectively-invariant characterization of conformal sructures corresponding 
to linearizable equations: the linearizing transformation is any projective transformation which 
sends H to the hyperplane at infinity (see Sect. 7.2). 

The tensor a^k provides a simple characterization of Lagrangian equations: 

Proposition 2. Equation ji^| ) is Lagrangian iff the corresponding conformal structure satisfies 
the following two properties: 

(1) the tensor is totally symmetric; in fact, since a^k is manifestly symmetric in the first 
two indices, it is sufficient to require a^k = a^j, 

(2) the covector Sj is a gradient: djSi = diSj; 

these conditions are obtained by weakening the corresponding conditions of Proposition 1 (recall 
that any linearizable equation is automatically Lagrangian). 
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Proof: 

To show that a given equation is Lagrangian, one has to find an integrating factor r such that 
the matrix rfyj is the Hessian matrix of a function, equivalently, dk{rfij) = dj(rfik), which 
gives 

fij &kfij — fik &j fik- 
T T 

Contracting this expression with f lJ one gets 

But 

J ^ = ^(d j fik-d k f ij ). (46) 
Substituting this back into the previous equation one obtains the relation 

jpq jpq 
dkfij fij~ 7 \®qfpk &kfpq) = 9j fik fik~ 7 {^qfpj 9j fpq) ; 

which is identical with = a^j. Finally, the right hand side of (I46D must be a gradient. Since 
the expression f^d^fij is automatically a gradient by virtue of the identity 

trF^dkF = d fc lndetF, 

one has to require that f l3 djfik is a gradient. This, however, is equivalent to the requirement 
that the covector Sj must be a gradient. This finishes the proof. 

Finally, we have the following 

Proposition 3. Conformal structure fijdp l dp ) in P n possesses a quadratic complex of null 
lines iff the symmetrized tensor a^ vanishes: 

afofc) = o. 

Indeed, the condition auj^\ = is identical to 



7.3 Tensorial formulation of the integrability conditions 

Let us begin with a general differential-geometric digression. Given a metric fij and a connection 
V with Christoffel symbols T l - k on a n-dimensional manifold, let us introduce the following 
objects: 

— covector s Sk and C&: 

°k = (n+2K»-i) (( n + 3 )^kfij ~ 2(n + l)Vjf ik j ; 

— tensor a^j.: 

O-ijk = ^kfij ipk ~l~ ^Skjfij &ifkj Sjfki] 

- tensor a ijk i: 

&ijki = V;Ojjfc — (q + 3si)aijk — Sjdy& — 9 jam. — SfcOyj. 

The importance of these objects is explained by their transformation properties: suppose that 
the metric fij and the connection V are allowed to vary within their conformal and projective 
classes, respectively, that is, 

fij ><Pfij, Fjk^Pjk + M + tA- 
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One can readily verify the transformation properties 

/ _l d M 

Sk —> Sk — Wki c k — > Cfc H — , dijk — > ^ijfcj OyfcJ — > ipdijkU 

which, in particular, imply that the new Christoffel symbols, 

r)k = r} fc + s fc <5j + ai 4, 

give rise to a well-defined affine connection V which depends neither on the choice of T l j k in its 
projective class, nor on the conformal factor ip. The expressions for Ojjfc nd ajj/y compactify to 

O'ijk — ^kfij Ckfiji Q>ijkl — V/Ojj'fc QOjjfc- 

In our context, n = 3, /y is a conformal structure in P 3 , and the projective structure is associated 
with the projective class of a flat connection: T l - k = (indeed, only projective transformations 

preserve the projective class of a flat connection). Thus, = d/dp k , and the connection V is 
given by 

T i jk = a k 8 i j + 8 j S i h ; 

notice that, although V is manifestly projectively flat, is does not need to be flat (that is, have 
zero curvatute tensor) in general. The tensors aijk, a-ijkl an d the affine connection V constitute 
a complete set of projective invariants of the conformal structure fijdp l dp ) . For n = 3 the 
integrability conditions (I13p can be formulated as follows (we begin with the Lagrangian case, 
which is computationally simpler): 

Integrability conditions in the Lagrangian case 

In the Lagrangian case the tensor a^k is totally symmetric, and the integrability conditions take 
the form 

djSi SiSj = '^Q'tmfiQ'rsvfiqfjpf^ ^ ^ "^j^Hpr^'jqsf^f ; (^7) 

and 

9 



aijki =j^Sym a ijp akiqf pq 



9 3 3 

SytTl | -^^fsl a ipqO'jtmfrk ~\~ ~^ltm^kpqf rif sj ~^^ a tmfj,0-rsufiqfjpfklf ) ^ & 



(48) 



respectively. Here Sym denotes a complete symmetrization in k, 

Sym T ijk = — T u(i)a{j)a{k) , 

and is the totally antisymmetric tensor dual to the volume form of the metric fy, that is, 
£ 123 _ \ jyfp^ £ 2i3 _ —\j^fp^ e tc, F = det/y. This provides yet another form of the integrability 
conditions in the Lagrangian case, compare with |20j . Recall that linearizable equations are 
characterized by the relations a^k = 0, djSi — SiSj = (Prop. 2 of Sect. 7.2), which clearly 
annihilate both of the conditions (|47p . (I48p . This is in agreement with the obvious fact that any 
linearizable equation is automatically integrable. 

Integrability conditions in the general case 

In the general case the tensor Oyfc is no longer symmetric (only in the first two indices), and the 
integrability conditions become considerably more complicated. Thus, the analogue of Eq. (|47f) 
takes the form 

1 r pq rs 

2q (pQ'iprQ'jqs ^Oyprittqsj + Q'iprQ'qsj Q>pri®'jqs) f ^ f j 



30 



one can show that the right hand side of (|49p is symmetric with respect to i and j, so that 
djSi = diSj. This means that, for integrable equations, the covector Sj must be a gradient. In 
this case the left hand side of Eq. (|49p can be represented in the form djSi — SiSj = ^rjRij 
where Rij is the Ricci tensor on V. The analogue of Eq. (I48p takes the form 

aijki = --^.Sym((8aki q + 4Aa qkl - 70a q ik)a>ijp + (64a w? - 82a qk i + 30a q i k )a pi j - 12a pk ja q u)f pq 



here Sym denotes symmetrization with respect to i and j. Both conditions (|49p . ()50p simplify 
to (|47p . (|48p under the Lagrangian assumption. These conditions provide a straightforward 
computer test of the integrability for equations from the class ([T]). 

7.4 Integrable equations and conformal structures possessing conjugate null 
coordinate systems 

Our first result is the following 

Theorem 5 The conformal structure fijdp l dpP corresponding to any integrable three-dimensional 
equation ([!]) is conformally flat. 

The proof is a straightforward calculation of the corresponding Cotton tensor, based on the 
integrability conditions derived in Sect. 2. Recall that for three-dimensional Lagrangian systems 
this result was established earlier in |20] . We emphasize that the transformation which brings 
the metric to a constant coefficient form is not necessarily projective (it does become projective 
for linearizable systems only). Thus, the theory of integrable equations of the form ([1]) has two 
'flat' counterparts: the first one is a flat projective structure provided by the projective space 
P n with coordinates p 1 , ...,p n , and the standard projective action of SL(n + 1, R). The second is 
the flat conformal structure fijdp l dp J . Although, viewed separately, both structures are trivial, 
this is no longer true when they are imposed simultaneously: their 'flat coordinate systems' do 
not coincide in general. 

Our next goal is to provide a differential-geometric characterization of hydrodynamic reduc- 
tions. Although our discussion will be restricted to the dimension three, all conclusions hold 
in any dimension. We will follow the notation of Sect. 2. Let p = (p l ,p 2 ,p^) = (a,b,c) be 
functions of the Riemann invariants R % (for our purposes it will be sufficient to consider one-, 
two- and three-component reductions only). By virtue of (j53|) . the derivative of p with respect 
to R l is given by 



The dispersion relation (|54p implies that dip a null vector of the conformal structure fijdp t dp J . 
Thus, the Riemann invariants R 1 provide a net of null curves on the corresponding submanifold 
p(R). Furthermore, the relation 




^Ojfirnt^urs)fipfj q fklf^ ~i~ 



{Q4:Cik q pfsj l^afcpfjfgj + 10cip q j f s k 2Qcijqpf sk )a,it m f r i-\- 

(102o,jgpO^ m j f>\a,p q iO,i m j A8cii q pCLj m t) frkfsl~\~ 




(50) 



d iP = (i,/AA l ) dm. 



didjp G span{9jp, djp} 
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which readily follow from ([7|), implies that this net is conjugate. Thus, we have the following 
geometric picture: 

One-component reductions correspond to null curves of the associated conformal structure. 
Two- component reductions are in one-to-one correspondence with surfaces which carry a 
conjugate net of null curves. Notice that we are in the realm of two different geometries, 
namely, conformal geometry (responsible for the property of being null), and projective geometry 
(responsible for the property of being conjugate). 

Three-component reductions correspond to three-conjugate null coordinate systems in P 3 . 
Since the existence of 'sufficiently many' three-component reductions is a necessary and sufficient 
condition for the integrability, the problem of the classification of three-dimensional integrable 
equations of the form ([T]) can be reformulated geometrically as follows: classify conformal struc- 
tures in the projective space P 3 which possess infinitely many three- conjugate null coordinate 
systems parametrized by three arbitrary functions of one variable. 

It is a truly remarkable fact that the moduli space of such structures is only 20-dimensional! 

8 Appendix: proof of Theorem 1 

Here we provide further details of the proof which was only sketched in Sect. 2. Our starting 
point is the quasilinear representation of Eq. ([1]), 

a y = b x , a t = c x , b t = c y , f n a x + f 2 2b y + + 2f V2 a y + 2/i 3 o t + 2/ 23 6 t = 0. (51) 

Following the method of hydrodynamic reductions, we seek multi-phase solutions in the form 

a = a(R\...,R N ), b = b{R l , . . . , R N ), c = c{R l , . . . , R N ), (52) 

where the phases R l (x,y,t), R N (x,y,t) are arbitrary solutions of Eqs. (J6j). Substituting 
the ansatz (|52p into (|5ip one obtains the equations 

dib = f^dia, die = A^a, (53) 

along with the dispersion relation 

I>(AV) = hi + / 22 (/i 1 ) 2 + / 33 (V) 2 + 2/ 12/ / + 2/ 13 V + 2/ 23 //y = 0. (54) 

Hereafter, we assume the conic (I54p to be irreducible. This condition is equivalent to the non- 
vanishing of the determinant of the coefficient matrix 

(fit fl2 fl3 \ 
fu hi /23 J , det F jk 0. 
fl3 f%3 /33 / 

The consistency conditions for Eqs. (|53p imply 

d A* d 

didja = —7^ — r^dia + ~ ■ - — -rdna. (55) 

Differentiating the dispersion relation (|54p with respect to R 3 , j 7^ i, and keeping in mind Eqs. 
(f5"3"P and ([7]), one obtains explicit expressions for djX 1 and djfi 1 in the form 

OjX 1 = (A* - X j )B ij d j a, djfj = (p? - n j )Bijd ja , (56) 
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where B^ are rational expressions in A* , X J , ^ , whose coefficients depend on fij(a,b,c) and 
first order derivatives thereof. Explicitly, one has 

B = % = i ^ . 

lJ Dij 2 f n + /aa/zV + + / 12 (^ + /*?) + / 13 (A< + Ai) + / 23 (^ + ^A*) ' 

notice that, modulo the dispersion relation (f54"|) . the denominator Dij equals AD ^ X% \ X3 , M 1^ 
The numerator iV« is a polynomial expression of the form 

Ntj = /33, 3 A i2 A^ + / 3 3,2^'A i2 + / 33 ,iA i2 + 2/23,3/^ V + 2/23,2^* + 2/23,1^ 
+/ 22 ,3^ 2 A J + /22,2/^V + / 22 ,i^ 2 + 2/i 3 , 3 A i A-'' + 2/ 13 , 2 ^'A i + 2/03,1 A* 

+2/12,3/^ + 2/i 2 , 2 /V + 2/i 2 ,i^ + /n, 3 A J ' + /11 V + /in; 

we adopt the convention that variables a, 6, c correspond to indices 1,2,3: thus, f\\ t \ = fii >a , 
etc. Taking into account Eqs. (fo^j) . Eqs. ([53]) assume the form 

didja = — {Bij + Bji)di,adja. (57) 

The compatibility conditions d k dj\ l = djd k X l , d k djfi l = djd k n l and d k djdia = djd k dia are 
equivalent to the equations 

dkBij = (BijBkj + BijBik — B k jB ik )d k a, (58) 



which must be satisfied identically by virtue of Eqs. (|53j) . (|o4j) . ([56]) and (|5Tj) . In order to obtain 
equations with 'simplest possible' coefficients at the second order derivatives of fij(a,b,c) we 
rewrite Eqs. (|58|) as 

d/ciVij = N{j (^-d k Dij + B kj d k a + B ik d k a^ - DijB k jB ik d k a. (59) 

The second order derivatives of fij(a,b,c) are present only in the l.h.s. term d k Nij. Further 
reduction of the complexity of the expression in the r.h.s. is achieved by representing 1/Dij in 
the form 

= Uij = [2(A% + /i 2 )(A j / 23 + /12) - /22(A l 'AV 3 3 + (V + A'")/is + / u ) 

+/ 22 (AV 23 + /12)// + / 22 (AV 23 + fu)v j + f&M/M* - A^') 2 detF), 

which holds identically modulo the dispersion relation (|54|) , and a subsequent substitution B s t = 
Ngt/Dgt = N st U s t. The denominators of the r.h.s. terms in Eqs. (j59j) cancel out, producing a 
polynomial in A 1 , A- 7 , A fc , with coefficients depending on the functions fij(a,b,c) and 

their derivatives up to second order. This was the most essential technical part of the calculation: 
the original expression (|59p has more than 1.000.000 terms with different denominators; after 
properly organized cancellations it reduces to a polynomial expression with less than 4500 terms. 
Using the dispersion relation (j54"|) and assuming, say, / 22 7^ (this can always be achieved by 
a linear change of the independent variables x,y,t), we simplify this polynomial by excluding 
the powers (/x l ) s , (^) s , (^ k ) s , s > 2, arriving at a polynomial of degree one in each of //, fJ, 
fi k , and degree two in A's. Equating similar coefficients of these polynomials in both sides of 
Eqs. (j59j) . we arrive at a set of 45 equations for the derivatives of the coefficients fij(a,b,c), 
which are linear in the second order derivatives. One can verify that only 30 of these equations 
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are linearly independent. Solving them, we get closed form expressions for the second order 
partial derivatives of the coefficients /n, /12, /13, /23, /33 in terms of the first order derivatives 
thereof, 30 equations altogether (without any loss of generality one can set = 1)) which can 
be represented in symbolic form (|13ft . 

d 2 fij = —R(fki,dfki); 

here R is quadratic in both and dfki- A straigtforward computation shows that this system 
is in involution: all compatibility conditions are satisfied identically. Since the values of the 
five functions /n, /12, /13, /23, /33, and first order derivatives thereof, are not restricted by any 
additional constraints, we obtain a 5+3-5 = 20-dimensional moduli space of integrable equations. 
This finishes the proof of Theorem 1. 
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